The Urge to Merge: When Cellular Service 
Providers Pool Capacity 

Sha Hua, Pei Liu and Shivendra S. Panwar* 



Abstract — As cellular networks are turning into a platform for 
ubiquitous data access, cellular operators are facing a severe data 
capacity crisis due to the exponential growth of traffic generated 
by mobile users. In this work, we investigate the benefits of sharing 
infrastructure and spectrum among two cellular operators. Specif- 
ically, we provide a multi-cell analytical model using stochastic 
geometry to identify the performance gain under different sharing 
strategies, which gives tractable and accurate results. To validate 
the performance using a realistic setting, we conduct extensive 
simulations for a multi-cell OFDMA system using real base station 
locations. Both analytical and simulation results show that even 
a simple cooperation strategy between two similar operators, 
where they share spectrum and base stations, roughly quadruples 
capacity as compared to the capacity of a single operator. This 
is equivalent to doubling the capacity per customer, providing a 
strong incentive for operators to cooperate, if not actually merge. 

I. Introduction 

The recent commercialization and popularization of 3G 
networks have significantly enhanced the mobile users' capa- 
bility of making use of ubiquitous data access. However, as 
users "mobilize" all of their communication activities, from 
streaming video to cloud computing tasks, cellular operators 
are facing a severe data capacity crisis. A recent study of 
Cisco predicts that mobile data traffic will skyrocket by a 
factor between 25X to 50X by 2015 [1]. Such explosive data 
traffic growth will overload the cellular infrastructure and result 
in either poor or expensive service to the subscribers. To 
address this challenge, the most straightforward solutions, such 
as adding cells and increasing spectrum have become either 
expensive or inefficient options. 

Cooperation among operators is yet another solution which 
has not drawn much attention from researchers. As cellular 
operators are devoting most of efforts in expanding their 
respective networks due to the urgency of the data capacity 
crisis and market competition, it has become apparent that there 
exist huge variations in spectrum usage, channel quality and 
coverage in different operators' networks [2|. Such diversity 
generates plenty of cooperation opportunities, which can be 
exploited to improve their network performance, spectrum 
efficiency and user experience. 

One motivating example is shown in Figure [T] Two operators, 
Operator 1 and Operator 2 (OPi and OP2 for short), provide 
cellular service in one area. We assume their base stations (BSs) 
are deployed following a hexagonal layout (solid and dotted 
lines, respectively) and are interweaved. Without cooperation, 
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Fig. 1. Motivating scenario of cooperation between cellular operators. 

the two operators will serve their own users solely using their 
own spectrum bands. The edge users, such as those at points 
B, D and F in the left cell of OPi, will experience bad channel 
condition and strong inter-cell interference. However, if the two 
operators cooperate by allowing their users to flexibly connect 
to any BS that is closest to them, users of OPi at B, D and F 
will be served by the BSs of OP2 and have excellent channels 
to the BSs. Generally, the users of OPi in the area AABC, 
AECD and AEFA will enjoy performance gains; these are 
mostly edge users with low data rates. A similar effect happens 
to the users of OP2 as well, such as those in the area AAOC. 
As we can see, simple cooperation improves the capacity of 
both operators and results in a win- win situation. Depending 
on the levels of synergy, the operators can share each other's 
infrastructure, spectrum, or even merge into a single operator. 
Interestingly, in AT&T's proposed acquisition of T-Mobile, one 
of the reasons cited is the resulting increase in capacity 0. 

In this paper, we focus on studying the benefits brought by 
different levels of cooperation among cellular operators. We 
are interested in answering the following questions: In what 
ways can the operators cooperate? How much performance 
improvement can their cooperation bring? Only a few recent 
papers address this issue [4| 1 5 1 1 6 1. In [4], Chandramani et al. 
optimize the operators' aggregate payoffs given BS locations 
and user rates, and use game theory to discuss how to share the 
profits. Supratim et al. in develop a user choice algorithm 
with network information provided by the operators. Peng et al. 
in (6) focus on spectrum-based cooperation, and form a group 
bargaining model based on the demand on each BS. 

However, all the above papers require cross-layer optimiza- 
tion based on the instantaneous network status. They do not 
characterize the long-term network capacity improvement for 
a large-scale multi-cell deployment. Our research fills this gap 
and provides a tractable and accurate model for quantifying 
the gains under different cooperation strategies. To validate 
our model, we have collected a set of real BS location data 
and provide performance results using real- world layouts under 



an Orthogonal Frequency Division Multiple Access (OFDM A) 
system. This aspect also has not yet been examined before. We 
expect our paper to give an insight into the benefit of cooper- 
ation and guide operators considering a range of cooperation 
options. Specifically the contribution of the paper is two-fold: 

• We present an analytical model for identifying the perfor- 
mance gain from cooperation among cellular operators. It 
provides a tractable and reasonably accurate model for av- 
erage user rate/throughput under a multi-cell environment. 

• We perform extensive simulations using real BS locations 
and OFDM A resource allocation algorithms. It validates 
the advantage of operator cooperation in a realistic setting. 

The remainder of the paper is organized as follows. Section 
II presents two modes of cooperation and the analytical model. 
In Section III, we provide an OFDMA resource allocation 
algorithm. The performance evaluation is presented in Section 
IV. Finally, the paper is concluded in Section V. 

II. Analytical Modeling 

In this section, we derive an analytical model to evaluate 
the network performance when the cellular operators cooperate. 
We consider two strategies, with different levels of cooperation, 
that we envisage the cellular operators may adopt: 

1) FLEXROAM (short for "flexible roaming"): the cellular 
operators allow their users to freely connect to a BS of either 
operator that provides the best signal strength. It will require 
an update in the signaling protocols to facilitate this. 

2) MERGER: in addition to FLEXROAM, all the BSs of an 
operator can reuse the spectrum of its cooperating operators. 
This may be the result of a full merger of the two operators, 
but does not preclude a mutual business agreement short of a 
full merger. 

Note that FLEXROAM and MERGER can also be offered to 
users by a mobile virtual network operator, or MVNO, which 
can purchase wholesale use of the infrastructure from the two 
operators. We further use NOCOOP to refer to the scheme that 
the two operators do not cooperate. 

Cellular networks are traditionally modeled by assuming the 
BSs are placed following a hexagonal layout Q. However, 
these models have long suffered from being both intractable and 
highly idealized. Recently, a general model based on stochastic 
geometry was proposed in (8), which provides tractable ways 
to evaluate network performance considering inter-cell interfer- 
ence and fading. Our analytical modeling follows the methods 
in 0, and provides tractable results of the performance under 
the various modes of operators' cooperation. 

We will mainly focus on the user's average ergodic rate 
and average throughput, which are important metrics for 
benchmarking the cellular system. We assume there are two 
operators, OPi and OP2, coexisting in the network. Our models 
can be easily extended to include more than two operators, and 
to derive other performance metrics, such as the distribution of 
user's Signal-to-Interference-plus-Noise Ratio (SINR). 

We assume the locations of BSs of OPi and OP2 follow 
independent Poisson point processes $1 and $2 with densities 
Ai and A2, respectively. The model drawn from such a random 



deployment is shown to be about as accurate as the standard 
grid model, compared to an actual cellular network [8 ]. Further 
we use Wi to denote channel bandwidths used by the BSs of 
OP^ i G {1,2}. We consider the downlink performance of the 
system. Without losing generality, we assume a typical user is 
located at the origin. If it is connected to the BS bo of OP^, its 
SINR can be expressed as: 

Pth rr a 
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where P t is the fixed transmission power of all the BSs, A^o is 
the noise power density. The distance between the user and BS 
bo is r5 , and the channel fading is h^. The user's distances 
to other interfering BSs b in &i are r 5 and the corresponding 
channel fadings are a is the path loss exponent. We assume 
all the fadings are Rayleigh fading with mean 1. 

A. Average User Ergodic Rate 

Following the proofs in [8], when the two operators do not 
cooperate and run their service as is, the average ergodic rate 
in the downlink for a typical user of OP^'s is given by: 
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Eq. ([2]) and ^ can be computed numerically and are used 
for performance comparisons in this paper. Next we examine 
the performance under cooperation strategies. When OPi and 
OP2 share infrastructure following the FLEXROAM strategy, 
we have the following theorem: 

Theorem 1: Under FLEXROAM strategy, for a typical user, 
the average ergodic rate is given by: 
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and F(Xi) is given by Eq. (3). W = {Wi,W 2 } and A = 
{Ai, A2} are the vectors for bandwidths and BS densities. 

Proof: We use r to denote the distance between the 
user and its associated BS. Under FLEXROAM, the user will 
connect to the closest BS from either of the operators that gives 
the best average signal strength. The union of all the BSs is still 
a Poisson point process & with density Ai + A2. The cumulative 
density function (CDF) of the distance r is then: 

F r (R) = Pr[r < R] = 1 - Pr[r > R] 

= 1 — Pr[no BS from any operator closer than R] 

= 1 _ e -(Ax + A 2 )^ 2 _ (6) 



Thus the probability density function (PDF) of r is: 

2 2tt(Ai + A 2 )r. 
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divided by the number of users in the cell. Under NOCOOP 
scheme, it is: 
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Moreover, we have 



Rflexroam(W, A) = / / r (r)E M> *(ln(l + SINR))dr. 
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We further use Xi = 1 to represent the event that the user is 
associated to the BS of OP;. Note Xi G {0,1}, Vi = 1,2 and 
X\ + X2 = 1. According to the rule of total probability, 
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wherein, SINR^ is given by Eq. (1). Using Di to represent the 
user's distance to the closest BS from OP^, in a manner similar 
to Eq. (7), the CDF of A is f Di (r) = e~ x ^2nXr. We have: 

Pr[Xx = 1] = Pr[L>i < D 2 ] 

= e- X27rr h7r\ 2 r 2 e- Xl7rr h7r\ 1 r 1 dr 2 dr 1 
Jo J ri 



Ai/(Ai+A 2 ). 



(10) 



Similarly, we have Pr[X 2 = 1] = A2/(Ai + X 2 ). Plugging 



Eq. ([9]) and ( |T0| ) into ([8]), and following the rest of the steps as 
in Appendix B of |8 ] completes the proof. ■ 

We next discuss the average ergodic rates when the 
MERGER strategy is employed. Under this scheme, each BS 
reuses the spectrum owned by the other operators, e.g., they 
can operate on the whole spectrum W\ + W 2 . As a result, the 
BSs from both operators interfere with each other. We then 
have the following corollaries from Eq. (2) and (3): 

Corollary 1: Under the MERGER strategy, for a typical 
user, the average ergodic rate is: 

Rmerger(W, A) = RNOCOOpiW! + W 2 , Ai + A 2 ). (11) 

Proof: This is same as the single operator case with BSs 
of density Ai + X 2 operate on bandwidth W\ + W 2 . We then 
have the corollary following Eq. (2) and (3). ■ 

B. Average User Throughput 

We will now determine the average throughput for a typical 
user under FLEXROAM and MERGER. We use rji to denote 
the subscriber density of OP^, i G {1,2}. Generally, subscriber 
density is far larger than the macrocell BS density, thus we 
assume r\i > Aj,Vz,j G {1,2}. 

According to the Law of Large Numbers, in NOCOOP, if the 
two operators serve their users independently, the number of 
users in each cell of OP^ will be approximately ^| = 
where S is the total area under consideration. We further 
assume that the scheduling decisions on the BSs ensure pro- 
portional fairness, which is a common design objective of all 
current and next generation cellular systems. With proportional 
fairness, the average throughput of a user from OP^ is its rate 
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Moreover, under FLEXROAM and MERGER strategy, the 

number of users in each cell converges to • Therefore, 

A1+A2 

the average throughput for a typical user is: 

Th FLEX ROAM(W, A, 77) = RflexroamCW, A)^±^ 

= -^-(AxiJ^Ai, A 2 ) + A 2 i? 2 (Ai, A 2 )). (13) 

where 77 = {771, ry 2 }. Finally, under the MERGER strategy, the 
average throughput for a user is: 
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Compared to the traditional hexagonal model that requires 
complex system-level Monte-Carlo simulations, our analytical 
model can be computed using a numerical tool in a short time. 

III. An OFDMA-based Multi-Cell System 

In the previous section, we described an analytical model for 
evaluating the network performance under different cooperation 
strategies of the operators. However, the model simplifies the 
real system in three aspects: 1) it assumes Poisson random BS 
deployment; 2) it uses the entire spectrum without considering 
subchannelization; 3) it assumes perfect resource allocation 
following the proportional fairness criterion. 

Due to the above constraints, it is still good to validate 
what the network performance will be in a practical multi-cell 
system under the various cooperation strategies. To evaluate the 
benefits of operator cooperation in realistic system, we plan to 
combine real BS location data and practical user- subchannel 
scheduling algorithms. As OFDM A is widely adopted by all 
the next-generation cellular standards (WiMAX, LTE) to divide 
the spectrum into subchannels, in this section, we present a low 
complexity OFDMA multi-cell resource allocation algorithm. 
Results with real BS locations are given in the next section. 
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Fig. 2. A typical OFDMA frame in a temporal-frequency domain. 

A typical OFDMA system partitions the radio resource in 
both frequency and time domains as shown in Figure [2] The 
total bandwidth is divided into C orthogonal subchannels, and 
one frame is divided into T time slots. A subchannel and 
time slot combination, termed a tile, is the minimum unit for 
resource allocation. In one BS, a tile can only be assigned to 
one user to avoid intra-cell interference. For each BS b, we 
use Mb to denote the set of users associated with it. Thus 
there are totally \Mb\ users in this cell. For a specific BS, 



Mb may vary depending on the different cooperation strategies 
the operators use. Specifically, in NOCOOP, Mb are the users 
from the same operator of b that are connected to b, whereas 
in FLEXROAM and MERGER, M b include the users from 
the other operators that treat b as the best BS. In addition, 
the number of subchannels C may also vary according to the 
cooperation strategy. It will increase in MERGER as each BS 
reuses the spectrum owned by the other operators. 

In the multi-cell scenario, neighboring cells may reuse 
the same tile depending on the inter-cell interference. As a 
result, the optimal multi-cell resource allocation problem is 
combinatoric and generally regarded as an NP-hard problem. 
Only heuristic methods exist so far l9llfT0ll[TT1l . Next, we 
present a centralized greedy algorithm to achieve a sub-optimal 
solution with low computational complexity. It assumes the 
users' channel state information is available at a Radio Network 
Controller (RNC), which coordinates all the BSs. Note that 
since our purpose is to quantify the gain brought by the operator 
cooperation, designing more practical algorithms is beyond our 
scope. Our algorithm is generalized from the ones in l9lll[TQl . 
It runs for each frame and aims to maximize the total utilities 
of all the users throughout all the cells in a frame. The utility 
of user i is defined as U(Ri) = log(Ri), with Ri being its 
data rate in the current frame. Using log function as the utility 
has been proved to be equivalent to finding a proportional fair 
solution lfT2t Assuming there are B BSs and N users, the 
complexity of the algorithm is 0(TCBN). 

Algorithm 1: Multi-Cell Resource Allocation Algorithm 

1 Rkxt — 0; YbxCxf — 0; X^xCxt = 0; 

2 for t = 1 : T do 
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for c : 



1 : C do 

- set of all the BSs; 
while ft / do 

6* ^— BS with the least number of assigned 
channels in Q; 

n <- ft 

Y' <- set Y b *, c , t = 1 in Y; 
foreach m G Mb* do 

L Atf m = Ug -Ul os 

m* <(— argmax m AUm\ 
if AUm* > then 

X i— set X m *,c,t = 1 in X; 
R <- \/i,Xi, c , t = 1, set 
Ri,t = Ri,t - n, c ,t(Y / ) + rv,t(Y) in R; 
Y <- Y'; 



As shown in Algorithm 1 , Line 1 initializes the matrices that 
contain the allocation results. Each entry R i t in R records 
the aggregate data rate of user i up to the slot t. X and Y 
are boolean matrices. The entry Yb, c ,t = 1 in Y if BS b is 
transmitting on subchannel c in slot t. The entry X^ Ci t = 1 
in X indicates that the user i is receiving on subchannel c in 
slot t. The subchannels are allocated sequentially (Line 3), and 
for each one, we check each BS for whether it should transmit 
on this subchannel in the current slot. To maintain fairness, 
the BSs with fewer number of assigned subchannels are given 



higher priority (Line 6-7). Our algorithm is marginal utility 
driven. For the current BS 6*, Line 9-13 picks the user with 
the maximal marginal utility AU. For each user m associated 
with 6*, AU m is defined as the difference between the utility 
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Y' is the matrix assuming BS 6* is transmitting on subchannel 
c in slot t. The function r m>c> t(Y) is the rate of user m on sub- 
channel c in slot t based on the allocation matrix Y. We assume 
the transmit power is equally split across the subchannels on 
each BS. Considering the inter-cell interference, r m ^ c ^(Y) is 
defined as: 



c ,t(Y) = W c \og 2 (l- 
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where Wc = W/C is the bandwidth of a subchannel, bo is 
the BS user m associates to. h\ m is the subchannel fading 
between BS b and user m on subchannel c, and db y m is the 
distance between BS b and user m. Finally, Line 12-15 update 
the matrices if the maximal marginal utility is positive. 

IV. Performance Evaluation 

In this section, we assess the performance of our proposed 
cooperation strategies through two sets of simulations. In the 
first set, we numerically compute the expressions given by 
the analytical model described in Sec. II. To further show the 
results under practical settings, in the second set, we conduct 
simulations for the OFDMA system with the multi-cell resource 
allocation algorithm shown in Sec. III. 

We have obtained precise coordinates for BSs from two 
major operators over a large suburban area near Washington 
D.C.. The 20 x 20 km area we chose is shown in Figure [3] 
There are 16 BSs from one operator and 13 BSs from the other. 
We mark the BSs using red and blue to distinguish them. Their 
location information will be used for our simulation. 
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Fig. 3. A 20 x 20 km view of the current BS deployment by two major 
cellular operators in the Washington D.C. area. 
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Fig. 4. The impact of BS density. 
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Fig. 5. The impact of user density. 
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Fig. 6. The impact of bandwidth. 



A. Numerical Evaluations 

Here we numerically compute the results of our analytical 
model when different cooperation strategies are employed by 
the operators. We will focus on user's average throughput as it 
is a critical metric for evaluating cellular service. It is important 
to show how much improvement there is through operator 
cooperation, and how such an improvement varies when the 
two operators have different BS densities, user densities and 
band widths. We set the value of A i in Eq. (12)-(14) as the red 
operator's BS density, e.g., Ai = 16/400000000 = 4 * 10~ 8 . 
Ai will be fixed in all the scenarios while A2 is left as an 
adjustable parameter. According to IEEE 802.16m evaluation 
methodology document 1T3TL we set the transmission power of 
the BSs as Pt = 46 dBm, noise power density as 7V = —174 
dBm/Hz and the path loss exponent a = 3.76. 

1 ) The Impact of BS Density: In this scenario, we set both 
OPi and OP 2 's bandwidths to W\ = W 2 = 10 MHz according 
to fT3ll . And we let their respective user densities be ^ = 
^ = 100. As a result, their spectrum resource and average 
number of users per cell are identical. Figure [4] shows the 
results by changing the BS density ratio A2/A1. As A2 changes 
from 0.2 Ai to 2A X , the performance of OP 2 under NOCOOP 
only slightly increases, since A does not influence the signal- 
to-interference ratio and the system is interference limited. 
We can easily see cooperation brings significant benefits to 
both operators. When ^ = 1, the two operators have 193.3 
kb/s average user throughput without cooperation. However, if 
they cooperate, FLEXROAM and MERGER strategies improve 
the average user throughput to 281.0 kb/s and 387.4 kb/s 
respectively, which is equivalent to a 45.4% and 100.4% in- 
crease, respectively. As A2 keeps increasing, the performance of 
FLEXROAM will slightly drop since the opportunistic diversity 
is reduced. However, FLEXROAM and MERGER still achieve 
a large performance gain for both operators within a wide range 
of BS density ratios. 

2) The Impact of User Density: To explore the impact of 
user densities on the network performance of cooperation, 
we fix A 2 = Ai = 4 * 10~ 8 and W 2 = Wi = 10 MHz. 
Then we let 771 = 100Ai and vary 772 from O.2771 to 2771. 
This represents different user/BS ratio for OP2, e.g., when 
772 = O.2771 = 2OA2, on average there are only 20 users per 



cell. Figure [5] shows the results. As we can see, as 772 increases, 
under NOCOOP, the average user throughput of OP2 drops 
drastically when more users are sharing the radio resource. In 
contrast, under FLEXROAM and MERGER the performance 
degrades gradually since users of OP2 can still opportunistically 
connect to the BSs of OPi. Moreover, when OPi and OP2's 
user densities are different, they achieve different performance 
sains. For instance when ^ = 0.8, under NOCOOP, OPi 
and OP2 have 193.26 kb/s and 241.58 kb/s average user 
throughput. When FLEXROAM is adopted, the performance 
jumps to 312.21 kb/s, which gives 61.5% and 29.2% gains to 
OPi and OP2, respectively (marked in the circles in the figure). 
Generally, FLEXROAM results in a win- win situation when ^ 
ranges from around 0.5 to 2.0. MERGER has a larger win-win 
range, which is from 0.37 to beyond 2.0. When the difference 
between 7/1 and 7/2 is too large, cooperation may generate 
negative effects to one operator. For example, when ^ = 0.4, 
OP2 loses performance under FLEXROAM. However, since we 
just are studying a simple roaming policy, such negative effects 
are expected to be avoided by adopting more advanced policies 
and incentive mechanisms. 

3) The Impact of BS Bandwidth: In this scenario, we fix 
771 = 772 = 100 Ai = 100A 2 . We fix W 1 = 10 MHz, and 
change the ratio W2/W1. As we can see from Figure [6] 
the performance of FLEXROAM, MERGER and OP 2 under 
NOCOOP increases almost linearly as the total bandwidth 
increases. When ^ > 0.4, FLEXROAM provides better 
performance for both OPi and OP2. It enables more gain for 
the operator with less spectrum. MERGER achieves a larger 
gain for all the points in the figure. 

B. Performance with Real BS Locations 

To validate the numerical results in practice, we conducted 
the OFDMA-based simulations with real BS locations. We 
consider two cellular networks with BSs placed as in Figure [3] 
Mobile devices are deployed uniformly in the 20 km x 20 km 
area with a density of 100 users per cell for each operator. The 
system is assumed to be coordinated by the OFDMA resource 
allocation algorithm we proposed in Sec. III. Based on the IEEE 
802.16m evaluation methodology document |13], the main 
system parameters are summarized in Table [T| For the channel 
gain, we model the fast fading component as Rayleigh fading 



TABLE I 
OFDMA System Parameters 



Number of subchannels 


32 


Number of slots per frame 


60 


BS Transmit Power 


46 dBm 


Noise power spectrum density iVo 


-174 dBm 


Channel bandwidth 


10 MHz 



with mean 1. We also considered shadow fading and model it 
as a log-normal random variable with a standard deviation of 
8 dB. The overall effect of antenna gain, cable and penetration 
loss and noise figure is set to dB. Finally, the path loss at 
the distance d is modeled as L(d)dB = 17.39 + 3.761og 10 d. 

Depending on the different cooperation strategies, the user 
set of each BS may be different. Under NOCOOP and 
FLEXROAM, each operator still aims to optimize the total 
utility of all the users connected to its own BSs, while 
under MERGER, we optimize the global utility of the two 
operators. We ran each simulation for 30 consecutive frames 
and we averaged the results of 5 simulation runs. Further, we 
plug the density values Ai = 16/400000000 = 4 * 1(T 8 , 
A 2 = 13/400000000 = 3.25 * 1(T 8 , 771 = 100Ai, rj 2 = 100A 2 
into our previous analytical model to resolve the numerical 
results for comparison purpose. 
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Fig. 7. network performance under OFDMA-based system. 

Figure [7] shows the results. Under NOCOOP, the perfor- 
mance of OP2 is worse than OPi since it has fewer BSs and 
the resulting average signal is weaker. Overall, the average 
throughput of all the users in the system under NOCOOP 
is around 263.8 kb/s. Cooperation substantially improves the 
performance of both the operators the real setting. Specifically, 
the user's average throughput increases to 386.24 kb/s in 
FLEXROAM and 519.31 kb/s in MERGER. These give 34.4% 
and 80.75% improvement to OPi and 68.21% and 126.16% 
improvement to OP2 respectively. Moreover, our analytical 
model accurately captures this trend and provides very close 
improvement predictions in terms of percentage. OFDMA 
system generally offers better performance than the numerical 
results since it better exploits the user/frequency selectivity. 
Also, the analytical model is pessimistic since it assumes the 
random deployment of BSs, which is discussed in (8). 

To further scrutinize the gain for the users, we draw the CDF 
of the average throughput of the users in one simulation run. 
Figure [8] shows the results in the domain of 0-200 kb/s. For 
NOCOOP, we use the statistics from all the users in the system. 
The curves validates the improvement result from cooperation. 
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For example, the median value of user's average throughput 
increases by around 40% and 100% under FLEXROAM and 
MERGER compared to NOCOOP. 

V. Conclusion 

This work investigates the potential benefits of cooperation 
among cellular operators. Using stochastic geometry, we pro- 
vide a tractable analytical model to derive the average user 
rate and throughput under two cooperation strategies. With 
real base station locations, extensive simulations over multi-cell 
OFDMA system further validate the performance improvement. 
We show even a simple cooperation policy with modest changes 
to existing networks can achieve around 30%- 120% capacity 
gains per customer under typical conditions. 
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